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Abstract 

We consider the question of removing the ultraviolet cutoff in a 2D Quantum Field 
Theory with an interaction term which is non-renormalizable by power counting. This 
model arises as the first non-trivial correction beyond the Gaussian approximation of the 
so called Capillary Wave or Drumhead Model, and is rather important from a physical 
point of view since it correctly describes the finite size effects of two-dimensional interfaces. 
Despite the fact that the interaction is non-renormalizable, we prove that for a large class 
of regularization schemes the finite and divergent parts can be separated in a simple way. 
Furthermore, the finite part is independent of the choice of cutoff prescription used. 
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1 Introduction 



A model which is widely used to describe fluid interfaces is the so called Capillary Wave 
Model (CWM) [Q]. In the literature, this model is also known as Drumhead Model. It is 
derived from the assumption that the fluctuations of the interface are described by an effective 
Hamiltonian proportional to the variation of the surface area with respect to the equilibrium 
solution (which is the solution which minimizes the area). In the case of an interface bounded 
by a rectangle of size Li x L2 the Hamiltonian has the following form: 

n = a dXi dX2 [V^(X1,X2)]2, (1) 

Jo Jo 

where a denotes the interface tension, and the field (j){xi, X2) describes the interface displace- 
ment from the equilibrium position as a function of the two longitudinal coordinates xi and X2- 
We shall adopt in the following periodic boundary conditions, i.e. 0(xi,X2) = 4>{xi + Li,X2) 
and (j){xi,X2) = (j){xi,X2 + L2). It is possible to generalize our results to other boundary 
conditions. 

Due to its non-polynomial nature, until recently the CWM has been studied only in its 
quadratic approximation. In this limit the Hamiltonian (|l|) becomes equivalent to a massless 
2D Gaussian model, see Q for details. Recently, in the first correction beyond the 

Gaussian approximation was studied in detail. Despite the fact that this corresponds to 
a scalar field theory with a non-renormalizable interaction term we found strong evidence 
that the resulting contribution to the interface free energy is universal, i.e. independent 
of the regularization scheme. We gave some general arguments supporting this assertion 
and also tested two different regularization schemes, finding the same cutoff independent 
contribution. In this paper we want to further pursue this analysis, by studying a wider class 
of regularization schemes. 

Let us mention two more reasons to be interested in this model. First, it can be shown 
by duality, that the Hamiltonian (|l|) can be used to describe the physics of Wilson loops 
in the confining regime of Yang-Mills theories. In particular, the model describes Wilson 
loops beyond the roughening transition, hence in the region (which is the most interesting 
from a physical point of view) in which the continuum limit can be taken. In this context 
the Gaussian approximation of the CWM Hamiltonian was discussed for the first time by 
Liischer, Symanzik and Weisz in Second, eq. (|l|) coincides with the Nambu-Goto model 
for the bosonic string (in a special frame). The corresponding quantum theory is anomalous. 
Depending on the quantization method one finds either the breaking of rotational invariance 
or the appearance of interacting longitudinal modes. However, it can be shown that in 
the infrared limit, i.e. when Li and L2 become large, these anomalous features disappear and 
the string theory becomes a simple Conformal Field Theory, with a well defined quantum 
behavior. This limit exactly coincides with the Gaussian approximation mentioned above. 

This paper is organized as follows: In section 2 we discuss the expansion of the CWM 
at the first order beyond the Gaussian approximation, to be called 2-loop contribution in 
the following. We claim that in the ultraviolet limit a finite part of this contribution can be 
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extracted in a natural way and is independent of the choice of the regularization scheme. This 
claim is made precise in form of a statement at the end of section 2. In section 3 we explicitly 
compute the cutoff dependence of the 2-loop contribution for a continuous momentum space 
cutoff and for the lattice cutoff scheme, thus providing proofs of our statement for these two 
important regularization schemes. 

2 2-Loop Expansion of the CWM 

Our aim is to evaluate the partition function 

Z = y" L>0 exp I -a J^^' dxi j^^ dx2 + [V^{xi,X2)Y | . (2) 

We rescale 4>{xi,X2) = 4>{^i,^2) / ^/^^ where Xi = ^iLi. Then we expand in the adimensional 
variable {aLiL2)^^ ■ The Hamiltonian (|l|) becomes 

n M = aL^L2 + Hg [4>] - -^-r'Hp M + 0({aLiL2)-^) , (3) 



with 



'0 Jo 

The 'asymmetric gradient' is defined through 



W.M = k I'd^i ['d^2[Vum,C2)? (4) 

JO JO 

np[<P] = U'dCi ['d^2{[Vum,C2)?y . (5) 



V„ = ( A , ^x-V2_|_ ] = (v„,i, V„,2) with u = L2/L1 . (6) 



dii d^2 

The corresponding (formal) expansion of the partition function eq. @) is 

Z = const • exp(-aLiL2) • • (^1 + ^-1^ d^i £ ^^2 ^([V„</.(ei, 6)]')') + •••)• (7) 
Here, we have introduced the Gaussian measure 

{O) ^ ^J^D(t> e^v{-\ dii j\i2[yum.i2)?) 0{(p). (8) 
The first order (1-loop) contribution Zi to eq. (|2|) is 

Zi= I^D<Pexp(^-^l\ci£dC2[Vum,^2)?) ■ (9) 
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The subscript a indicates that the functional integral (which is divergent in the ultraviolet) 
is to be regularized by a suitable cutoff. Examples for cutoff schemes will be given below. Zi 
is a purely Gaussian integral and can be evaluated by using standard techniques. The result 
turns out to be a function of the adimensional ratio u alone: 



1 -2 
Zi {u) = —^r] (iu) /r] {i) 
Ju 



where r\ is the Dedekind eta function 



Ty(r) 



1/24 



q = exp(27rir) 



(10) 



(11) 



n=l 



This is a well known result in string theory and conformal field theory, and coincides^ with 
the partition function of a 2D conformal invariant free boson on a torus of modular parameter 
T = iu p. 

The calculation of the 2-loop contribution is less simple. One can show that 



(K<^(ei,e2)]')") = (-i)" E ••• E 



92 



exp 



k=0 



Here, we have defined, for /i = 1, 2, 



(12) 



(13) 



By using eq. ( |12[ ) for n = 2 and translational invariance on the torus, we find for the 2-loop 
expansion 



Z = const • exp(— (TL1L2) ■ Zi ■ il + 



R 



aLiL2 



+ 



R=- + 3G^ + 2G1G2 



with 



R can be rewritten as 

R=^ (2(Gi + G2f + (Gi - G2)' 
We now write R as a, sum of two parts, 

R = Roo + Ri , 

with 



Rn 



Rf 



(Gi + G2)2 , 



-(G1-G2 



(14) 
(15) 
(16) 

(17) 

(18) 
(19) 



""^The constant 77 (i) has been introduced to normalize Zi{l) = 1. 
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With these definitions at hand, we now formulate our 



STATEMENT : 

At least for the four regularization schemes 

Ai : (^-function regularization, 

A2 : point splitting regularization, 

A3 : continuous momentum space cutoff, 

A4 : lattice cutoff, 



the following assertion holds: 



The coefficient R of the 2-loop contribution to the CWM partition function can he naturally 
split as R = Roc + Rf, and 

(SI) Roo diverges when the cutoff parameter a is sent to zero. The singularity depends on the 
chosen regularization, but in all regularization schemes it is "shape independent" . For any 
finite value of the cutoff, Rqo is proportional to the area L1L2 of the rectangle, but does not 
depend on the ratio u = L2/L1. 



(S2) Rf remains finite when the cutoff is removed: 



Rf 



1 



2 + ( -uE2 (iu) - 1 



E2 denotes the first Eisenstein series: 

i?2(r) = l-24^- 



00 r) 

n q 



q = exp(27rir) . 



(20) 



(21) 



n=l 



The calculation of R by means of the (^-function regularization can be found in |^ , while the 
analogous calculation with the point-splitting procedure was presented in . We refer to the 
original papers for the details. Before we turn to the proofs for the two schemes A3 and A4, 
let us state two immediate consequences of part (5*1) of the statement. First, the fact that 
Roo depends on Li L2 alone (and not on the shape factor u) implies that ratios of partition 
functions of interfaces with different shapes but the same area stay finite when the cutoff is 
removed. Second, this same fact allows to 'absorb' the infinity in a 'renormalization' of the 
interface tension a (see also Q). 



3 Proofs 

In this section we present the proofs of the above statement for the cutoff schemes A3 and 
A4. We return to dimensionful variables in this section. 
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3.1 Momentum Space Cutoff 

A very general cutoff procedure for fields on a torus is to use the Fourier representation and 
then cut off large momenta. In such a scheme the regularized version of can be written 
as 

GM = -E§x(aV). (22) 

The momenta are summed over 

27r . 

Pix = -j—- integer , (23) 

and = Pi + P2- The cutoff function x is assumed to have the property x(0) = 1 and to 
drop sufficiently fast to zero when the argument goes to infinity. For technical reasons, we 
shall later assume that there exists a constant c < such that 

lx(p)l<e^^ (24) 

In order to remove the cutoff we have to send a to zero. Let us first study Gi + G2- 

Gi + G2 = l-Ex(«V), (25) 
p 

where the sum is over all momenta (including 0) . We introduce a Fourier representation of 
the cutoff function 

x(aV) = J d^xe'P-U^). (26) 

Here, 

Ux) = J ^e-^'-x(aV). (27) 



Now, 



Ex(aV) = / d'xU^) Y^e'P^ . (28) 



p p 
The sum over the momenta can be performed: 

J2 e'"'' = LiJ2 ^(^1 - niLi) L2 J2 ^i^2 - "2^2) . (29) 

p m n2 

Hence we have 

^ x{aV) = L1L2 ^ Xa{niLi,n2L2) . (30) 

p n\ n2 

The Fourier transform of the cutoff function is 

U^) = /(|f2^"^"^x(aV) 

= 2^1 dppx{p'')Mp\x\/a). (31) 
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Jo denotes the Bessel function of order 0. Plugging this into eq. (|30|), we get 



(32) 



In the hmit a ^ all contributions but the ni = n2 = term vanish. The argument is that 
Jo oscillates wildly in the relevant interval, with the consequence that for nonvanishing 
the integral vanishes faster than a?. Since Jo(0) = 1, we end up with 



Ex(aV) 



L1L2 



L1L2 
47ra2 



(33) 



This completes the proof of the part (51) of the statement. 

Let us now turn to the discussion of Gi — G2. In this case we use the Laplace transform 
of X- Foi' function with 

\fit)\<e^' (34) 



there is a representation 



For the cutoff function, 



c+ioo 



x(aV) 



1 
27d 



c+ioo 



(35) 



(36) 



From the assumption eq. (|24| ) we infer that we can choose c < 0. So the integration is over 
a contour where the real part of s is always negative. This allows for a representation of 
Gi — G2 as follows: 



Gi — G2 



2m 



dsx{s)F{s,a) , 



with 



Fis) 



^ Pi -Vl ^sa^p^ 



p. 



(37) 
(38) 



Writing out explicitly the momenta, we get 



F(s,a) = -E 



ni,n2 



Li 



2 exp 

1 I "-2 



s 47r^a^ 



n 



(39) 



Here, the prime on the summation symbol indicates that ni = ?i2 = is to be excluded. Let 
us now separate the two terms in the numerator, and let us call them Pi and P2. 



Pi 



E 



exp 



ni,n2 ~r 



(40) 
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and 



E 



ni,n2 H 



(41) 



r2 I r2 

In Pi we can safely take the limit a ^ in the 2-direction, because, as a sum over the index 
n2, Pi converges absolutely. We find 



E 



"l I "2 

"1,^2 + 



^ 2 2 '^l 



(42) 



Here we have introduced the symbol = which shall signal equality in the limit a ^ 0. Eq. (42) 
can be rewritten as follows: 



oo oo 



2 ^2 



U n 



^1^-2 E E ;^2^^«-p 

ni=ln2=— oo ^ 



s47r a 



2„2 ™1 



(43) 



The sum over n2 can be done explicitly by using the identity (see, e.g., text books on quantum 
statistics) 



E 



vr ' — ' + 6^ 

712 = — oo ^ 



1 + 



,27rfe 



1 



for 6 > . 



We get 



Pi= -2ttuI ^ ni 1 + 



ni=l 



exp 



s47r a 



2„2 



Let q = e^'^*'^ = e Then we have 



Pi = — 27ru ^ [ n exp 



71=1 



n 



^ 1 1 - a" 



exp 



n=l 



s 47r^a^ 



(44) 



(45) 



(46) 



The second term is not singular in the limit a ^ while the first term can be treated with 
the Euler-McLaurin summation formula: 



k 

E/N 

n=0 



fc 1 oo o 

/(x)dx + -[/(0) + /(A;)]+ y(-l)™[/(2™-i)(fc)-/(2— i)(0)l (47) 

m=l 



where the are the Bernoulli numbers. Taking f{n) = ne~*"^ we get /(2™-i)(oo) = and 



(2m-l) 



yielding 



/ 



n=0 



(0) 



1 d' 



'2m 



dx xe 



h' 



1 



-1 



x=0 



ml 



R 



m=l 



m!(2m)! 



(48) 



(49) 
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The last sum is an absolutely convergent series which is regular for 6^0. The integral is 
elementary: 



xe-"^ dx = —. (50) 

^0 



Taking b = —s in a /Li, we get 



oo 

n exp 



n=l 



.2 2 



n 



LtLo tt 

+ + - , (51) 



sAtto? 6 



where the omitted terms have higher powers of a. The finite part of this sum is exactly what 
is needed to reconstruct the Eisenstein function (once the limit a — > is taken in the second 
part of Pi) and we find: 

If we now study P2 we see that we only have to exchange Li and L2 and change the sign. In 
this way the divergent parts cancel out: 

F(s, a = 0) = -uE2{iu) - -u~^E2{iu'^) (53) 
6 6 

Using the identity 

1 fir 

E2{ — )=T^E2{T)-i— (54) 
r vr 

we arrive at 

F(s,a = 0) = (««)-!• (55) 

Note that when we performed the limit a — > 0, the function F became also independent of s. 
As a consequence, the Laplace integral eq. (|37| ) then becomes just a representation of x(0)) 
which is one, and part (5*2) of the statement is proven. 

3.2 Lattice Cutoff 



In this case is defined as: 



The momenta are summed over 

27 

al 



G^, = -y . (56) 



2tt 

Ppi = —r-n^,, = 0, . . . - 1 , L^ = al^, (57) 



and 

aYu = 4 sm\ap^) . (58) 
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To remove the lattice cutoff one has to send li and I2 to infinity, while keeping the ratio 

- = — =u (59) 

fixed. Gi can be written as 

h-i 

n,=ilnf^osin2(^)+sin2 (^) 

The sum over n2 can be done explicitly. To this end we employ the identity 
sinh2(x/2) 



h-i ( I2-1 sin2 f 
Gi = -T It 7 , r \ ■ (60) 

^-^ I ^-^ ■ 9 / 7rT7.i \ , ■9 / irnn \ I ^ ' 



sinh2(x/2) +sin2(u;/2) 



tanh(x/2) ^ ex.p{—x\n\ — iojn) , (61) 



which can, e.g., be derived by a series of elementary operations from formula 3.613 of ||8|. We 
make the identifications 

UJ = 27rn2//2, 
sinh2(x/2) = sin2(7rni//i) . (62) 



Using that 



h — l h — i 00 

exp(— iwn) = ^ exp(— i 27r n n2/^2) = ^2 ^ 6{n — l2k), (63) 

n2=0 n2=0 fc=— 00 



we find after summing over n, 



„f;„sin2(32!l)+sm2(2g2) l-exp(-;2i); 

This allows us to identify the convergent and divergent parts in the limit li,l2 00. We 
write 

Gi = Ci+Di, (65) 

where 

Di = -h T tanh(x/2) , (66) 

ni=l 

and 

Ci = -2 h E tanh(x/2) ^ ^^P^-/^;^) . (67) 

l-exp(-/2x) 
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In all these sums the index ni is hidden in x. Di can be written as 

h-i sm(^' 



"1=1 ^l + sin2(^) 

This sum can be evaluated by using the Euler-McLaurin formula eq. (|4^) and gives 

D, = -'-l^ + lu + 0{u/ll). 

In C2 we can safely take the limit li,l2 00. In this limit, we can write x = 27rni//i+0(l//f ). 
We can thus approximate tanh(x/2) ~ vrni/Zi, and 

exp(-/2 x) = exp(-27rnni) (l + O (uAi)) . (70) 

Let us define as usual q = e"^'^". Then 

00 n 

C^ = -4TTuy (71) 



Here we have taken into account that the two values ni and h — ni are identified with the 
same value of x thus giving a further factor of 2. 

Obviously, the corresponding results for G2 can be obtained just by letting u u~^. 
Putting together the various parts, we get 

Gi = Di + Ci^--^ + ^uE2{m), 

G2 = D2 + C2^--^ + ^u-'E2{iu-'). (72) 
Employing again eq. (^), we now easily verify that 

G1 + G2 ^ -(/1/2-I), 

G1-G2 ^ ^E2{iu)-1, (73) 
which proves the assertions of the statement. 



4 Conclusions 

It could seem rather surprising that we find a universal result, independent of the regular- 
ization scheme, for the finite part of R, since it stems from a non-renormalizable operator. 
We think that the reason, as it was already mentioned in Q and is that the result is 
completely fixed by symmetry requirements. In particular, the modular invariance of the 
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original interaction term (which is a direct consequence of the fact that the model, due to 
the periodic boundary conditions on <j), is defined on a torus) and its scaling behavior (which 
immediately follows from power counting) tells us that Gi — G2 must be a modular form of 
weight 2. The normalization coefficients in front of the modular form then follow from the 
study of the one-dimensional limit u — > 00, which was already discussed by Arvis in [p. 

The contribution that we have found does not depend on the regularization schemes that 
we have discussed because they all preserve the modular symmetry. We are thus led to ask if 
there is some physical reason to choose among all the possible schemes those which preserve 
this symmetry. We have no definite answer on this point but it is tempting to conjecture that 
it is a consequence of the fact that the Hamiltonians (Q) and (^) are respectively the first 
and the second order expansion of a bosonic string action, and that the modular invariance 
is actually preserved at the quantum level in the string model. 
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